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STABILITY OF FREE-CONVECTION BOUNDARY- LAYER FLOWS 
By Philip R. Nachtsheim 


SUMMARY 

The stability of free- convection boundary- layer flows is investigated by 
numerical integration of the disturbance differential equations. Stability cal- 
culations are carried out for Prandtl numbers of 0.733 (air) and 6.7 (water) with 
and without temperature fluctuations. Results presented for these four cases 
consist of eigenvalues (phase velocity, wave number, and Reynolds number), eigen- 
functions, and energy distribution curves for neutrally stable disturbances. 
Tabulations of the basic velocity and temperature profiles for a Prandtl number 
of 6.7 are also included. 

When temperature fluctuations are included, a mode of instability is found 
in which the primary source of energy for the disturbance motion arises from the 
interaction of buoyancy forces with velocity fluctuations. The present stability 
results are compared with available theoretical and experimental results. 


INTRODUCTION 

It is well known that boundary- layer flows are unstable under certain con- 
ditions in the sense that a small disturbance imposed on the basic flow can grow 
indefinitely in time. This instability is related to the shearing motion of the 
basic flow and existence of a mechanism for transferring energy from the basic 
flow to the disturbance motion via the Reynolds stress. It is also well known 
that a static fluid in a gravitational field in which a constant temperature 
gradient is maintained is unstable under certain conditions. This thermal in- 
stability arises when the gravity vector has a component parallel to the tempera- 
ture gradient. In this case a supply of energy for the disturbance flow comes 
from the potential energy of the fluid particles in the static configuration. 

Both of these stability problems have received extensive, but separate, treat- 
ments in the literature. 

The problem of the stability of free-convection boundary- layer flows is an 
interesting combination of the problems of boundary- layer stability and thermal 
stability. Since the free-convection flow is a shearing motion, the problem of 
its stability contains all the features of the boundary- layer stability problem. 
However, when a temperature disturbance is imposed on the free-convection 
boundary- layer flow, the essential feature of the thermal stability problem is 
introduced, and there is a possibility of transferring energy to the disturbance 
flow if there is a correlation between the disturbance velocity component along 



the plate and the disturbance temperature gradient along the plate. 

The determination of the stability characteristics of a given free- 
convection boundary- layer flow requires the determination of the eigenvalues 
occurring in a system of complex differential equations of sixth order. This 
system of equations differs from the well-known Orr-Sommerf eld fourth-order 
equation, which applies to ordinary boundary-layer flows in that the sixth-order 
system takes account of the interraction of the gravity force with density fluc- 
tuations. 

The effect of temperature fluctuations has been neglected in all previous 
investigations of the stability of free-convect ion boundary-layer flows. Refer- 
ence 1 showed that the neglect of temperature fluctuations is justified if the 
Reynolds number is sufficiently large; however, the Reynolds numbers at which 
finite disturbances have been observed experimentally in free-convect ion flows 
are not extremely large numbers. 

With the neglect of temperature fluctuations the disturbance differential 
equation for free-convect ion boundary- layer flows reduces to the Orr-Sommerf eld 
equation. The asymptotic techniques that were developed to solve the Orr- 
Sommerfeld equation were used in references 1 and 2 to solve the free -convection 
boundary-layer stability problem. 

In asymptotic methods, the critical layer plays a central role in the anal- 
ysis. The critical layer is located at that point in the boundary layer where 
the phase velocity of an assumed disturbance wave matches the velocity of the 
basic flow. For conventional velocity profiles, linearly independent solutions 
of the Orr-Sommerf eld equation are often obtained in the form of expansions about 
the critical layer. Since the velocity in the free-convection flow is zero at 
the plate and at large distances from the plate, each value of velocity is at- 
tained twice, once when the velocity is increasing and once when the velocity is 
decreasing; hence, for these profiles there will be two critical layers or none. 
In references 1 and 2, account is taken of the presence of two critical layers 
by obtaining solutions to the differential equation by means of expansions about 
each critical layer. However, minimum critical Reynolds numbers, that is, the 
Reynolds numbers below which all disturbances should be damped, were obtained 
that were greater than the Reynolds numbers at which finite disturbances were 
observed experimentally (refs. 1 and 2). 

In order to avoid the problem of expanding solutions about two critical 
layers, other investigators have used direct numerical methods to solve the Orr- 
Sommerfeld equation. A numerical method to determine the stability characteris- 
tics of a transverse velocity profile near the stagnation point of a sweptback 
wing was developed in references 3 and 4. This method consists essentially of 
using step-by-step integration to find two linearly independent solutions of the 
differential equation that satisfy the boundary conditions at the wall. A suit- 
able linear combination of these solutions is then sought by a trial-and-error 
method in order to match the known free-stream solution at the edge of the bound- 
ary layer. In reference 5, the method developed in reference 6 was applied to 
determine the stability characteristics of a free-convection velocity profile 
without temperature fluctuations. A minimum critical Reynolds number below the 
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Reynolds number at which, finite disturbances were observed experimentally was 
obtained in reference 5. However, none of these methods made use of any of the 
analytic properties of the Orr-Sommerf eld differential equation in order to 
simplify the numerical calculations. 

In the present numerical method, use is made of the property of the dis- 
turbance differential equations that solutions of the differential equations 
depend analytically on the parameters appearing in the equations. This method 
of solution, unlike the asymptotic methods, does not rely on the existence of a 
critical layer. 

The basic idea of the present method is to integrate the differential equa- 
tions by a step-by-step method for guessed values of the eigenvalues and start- 
ing values as if they were nonlinear equations. The boundary-value problem is 
treated as an initial value problem for which not all the proper starting values 
are known so as to satisfy conditions at the other boundary. Equations involv- 
ing the dependent variables are formulated and evaluated at the edge of boundary 
layer so that the zeros of these equations indicate that the numerical solution 
matches the known free-stream solution. The Newfcon-Raphson second-order process 
is used to obtain corrections to the eigenvalues and the starting values in order 
to satisfy boundary conditions. The iteration process is continued until the 
boundary conditions are satisfied. 

This general method outlined is described In reference 7. The purpose of 
this report is to apply the method to the solution of the free-convection 
boundary- layer stability differential equations and to show how the analytic 
properties of the differential equations can be used to simplify the application 
of the general method. 


DISTURBANCE EQUATIONS AND THEIR SOLUTIONS 
Disturbance Equations 

The disturbance equations for a parallel flow as given in reference 1 are 
as follows: 

cp rlM - 2a 2 cp M + a^cp + s' = laRe^F’ - c)(cp" - a 2 cp) - F ,n cpJ (l) 

s" - a^s = iaRePrJjF' - c)s - cpH^J (2) 

(Symbols are defined in appendix A. ) 

Equations similar to these were also derived in reference 2 by means of a 
different nondimens ionalizing procedure. A brief derivation of equations (l) 
and (2) is given in appendix B. The effect of the temperature fluctuations 
appears in the function s, which couples equations (l) and (2). If s’ is set 
equal to zero, equation (l) reduces to the well-known Orr-Sommerf eld equation. 
The boundary conditions for equations (l) and (2), in the case of an isothermal 
plate, are 
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'Q = 0 


00 


cp = cp r = S = 0 
cp, cp' , s -* 0 


(3) 

(4) 


The boundary condition that the temperature disturbance amplitude s van- 
ishes at the surface is a reasonable one for metallic plates immersed in a fluid 
such as air or water, where the plate materials are highly conductive compared 
to the fluid. 


Replacement of Boundary Conditions at Infinity 

The boundary conditions at infinity can be replaced by different conditions, 
which are to be satisfied at a finite distance from the plate. The new condi- 
tions are to be applied at the edge of the boundary layer at r\ = r| e , where the 
values of F’, F MT , and H 1 are nearly at their asymptotic values. For free- 
convection profiles F', F™, and H ! approach zero at the edge of the boundary 
layer; therefore, for rj > rj e the coefficients in equations (l) and (2) can be 
treated as constants, and equations (l) and (2) reduce to 

cp T ™ - 2a 2 cp n + aAp + s' = -iaRec(cp ,T - a 2 cp) (5) 

p 

s M - a s = -iaRecPrs (6) 

The general solution of the system of equations (5) and (6) is 
cp = c 1 exp (at) 4- c 2 exp (-at) + c^ exp (pt) 

+ c 4 exp (-pt) + c 5 exp (rt) + Cg ex P (-"Pi) ( 7 ) 


s = 


- c 5 ( r 


2 \ / 2 q 2 \ 

a )(r - P ) 


r 


exp (rt) + c 6 


T. a . ) (x ft , ) , exp (-rt) (8) 


T 


where Cp Cg, . . . , c^ are arbitrary constants and where p and y, respec- 
tively, are the roots with positive real parts of 


p 2 - 


a - iaRec 


(9) 


r 2 = a 2 - iaRecPr (10) 

The task remains of eliminating the arbitrary constants in equations (7) 
and (8) in order to obtain a set of linear and homogeneous relations in the de- 
pendent variables and their derivatives, which are to be satisfied at the edge 
of the boundary layer. The elimination of the arbitrary constants is carried 
out in appendix C, where it is shown that solutions of equations (l) and (2), 
which decay exponentially as the distance from the plate increases indefinitely, 
are characterized by three conditions that must be satisfied at the edge of the 
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boundary layer, rj = r) . These conditions are 

s' + ys = 0 (11) 

<P ,M - a-V + p(cp" - a 2 cp) 4- - s = 0 (12) 

cp'" + aq>" - p 2 (c p' + acp) 4 J — a - 0 (13) 

T i oc 


Solution of Eigenvalue Problem 

The approach to the eigenvalue problem for fixed Pr used herein is to fix 
a and Re and then to find values of c = c r + ic^ (eigenvalues) for which 
equations (l) and (2) have solutions (eigenf unctions) that satisfy the boundary 
conditions (eqs. (3) and (11) to (13)). 

The eigenvalue c and the corresponding eigenfunctions are obtained by 
treating equations (l) and (2) as if they were nonlinear equations. A trial 
solution is obtained by step-by-step numerical integration of the differential 
equations starting at rj = 0 for assumed starting values and an assumed value 
of c. Equations (11) to (13) are evaluated at the edge of the boundary layer. 

If the boundary conditions are not satisfied, the starting values and c are 
adjusted, and another trial is made. For the assumed values of the parameters, 
denote the values taken on by the expressions on the left side of equations (ll) 
to (13) by 

A.^ = s' + rs (14) 

Ag = cp"'- a 2 cp 4- 0(cp" - a 2 ) 4- - ■ n ■ s (15) 

A 3 = cp m+ acp" - |3 2 (<p' + «P) + - | - a (16) 

It is desired to find the eigenfunctions and the values of the parameters 
that 'will simultaneously cause Ap A^, and A^ to take on the value zero. In 
order to calculate Aq, Ag, and A3, starting values cp(0), cp’ (0), cp"(0), cp m (0), 
s(0), and s f (0) have to be assigned as well as a value for c. The boundary 
conditions at the plate are satisfied by taking cp(0) = cp f (0) = s(o) = 0. Since 
the differential equations (l) and (2) are linear and homogeneous and the bound- 
ary conditions are homogeneous, it is permissible to take cp fT (o) equal to some 

convenient fixed value. This value sets the magnitude of the linear oscillation. 
The remaining unknown starting values, cp fn (0) and s 1 (0), are denoted by a and 
b, respectively. 

Newton-Raphson method is used to obtain the zeros of equations (14) to (16). 
If the chosen values a, b, and c produce a solution cp(r);a,b, c) and 
s(r;;a,b, c) that approximately satisfies equations (14) to (16), that is, actually 


5 



leads to values of Ap Ag, and A^ close to zero., then a better approximation 
is obtained by starting with a + Aa, b + Ab, and c + Ac instead of a, b, 
and c, where Aa, Ab, and Ac are obtained as solutions of 


0 


0 


0 


— AjRa, b, c) + 
= Ag(a,b, c) + 
= A 3 (a,b, c) + 


Sa. 

5a 


Sa Sa 

x Aa + ,_Jl Ab + Ac 


5^ 




Sa Sa Sa^ 

Aa + — Ab + ^_£ Ac 

Sa Sb Sc 


5a~ 


Aa + 


5b“ 


Ab + 


Sa 3 

SF" 


Ac 


(17) 

(18) 
(19) 


The quantities k V A g, an( i and the partial derivatives evaluated at 

rj = rj e can be determined by integrating the differential equations (l) and (2) 
by a step-by-step method. The partial derivatives in equations (17) to (19) may 
be given either by a finite difference approximation or from equations obtained 
by partial differentiation of equations (l) and (2) vith respect to a, b, 
and c. (The coefficients of eqs. (l) and (2) are analytic functions of r] and 
depend analytically on the parameters a, Re, and c. Therefore, the solutions 
of eqs. (l) and (2) have the same analytical properties and have the required 
partial derivatives. ) In the former method four integrations of equations (l) 
and (2) are carried out, one vith a basic set a, b, and c, and three more, 
each vith appropriate small increments on a, b, and c. Partial differentia- 
tion gives the same information vith one integration. Hovever, in this second 
method, three extra sets of differential equations have to be integrated. 

The equations for the derivatives vith respect to a are as follows: 


cp" M - 2 cgV' + ctAp + s' = iaRe(~(F f - c)(cp" - a 2 cp ) - F”'cp 1 
a aa *- aa a_j 

s" - a 2 s = iaRePr Rf 1 - c)s - cp H'l 
a a L a a j 

with the initial conditions 

ri = 0 cp a = cp^ = <J>|[ = s a = s a = 0 = 1 

The equations for the derivatives vith respect to b are as follovs: 

cp”" - 2a 2 cp" + aV + S’ = iaRe Rf’ - c)(cp’ - a 2 cp ) - F'”cp] 
b bbb L bb J 

s^ - cc 2 s fe = iaRePr [(F' - c)s^ - cp^H^] 

with the initial conditions 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 
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TJ = 0 


<Pb = ^ = V* = < Pb" = S b = 0 


■b" 1 


(25) 


and the 
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tut 

c 


equations for the derivatives with respect to 
- 2a 2 cp ir + a^cp + s T = iaRe("(F r - c)(cp" - a 2 cp ) 

C C C L C C 


c are as follows: 

- F"'cp c - (cp" - a 2 cp)] 


j" - a 2 s = iaRePr f(F' - c )s - cp H' - si 
c c L c T c J 


with, the initial conditions 


^ - 0 <P C = % = *0 = ^c" = s c = S C = 0 


(26) 

(27) 


(28) 


In these equations the subscripts a, b, and c denote partial differentia- 
tion with respect to a = cp'" (0), b a s' (0), and c, respectively. For example, 
cp a = (d/da)cp(r];a, b, c) gives the rate of change of cp with respect to a with 
b, c, and tj held constant as a function of T|. The partial derivatives of the 
A' s appearing in equations (17) to (19) can now be evaluated. Note that the 
boundary conditions (eqs. (14) to (16)) involve c, the eigenvalue, through the 
intermediate variables 3 and y defined by equations (9) and (10). Utilizing 
equations (14) to (16) results in 


dA^ 

da a 


= s' + rs 

o 1 c 


dA. 

5V ’ 34 + Ts * 


s — = s' + r® 
oc c 1 c 


idRePr 1 

2 r s 


5 a 


= cp a " - cc 2 cp D + 3(<p" - a 2 cpj + 


r a 


r + 3 a 


dAg 

5b~ = <c - a \ + 0(<Pb - a \ ) + y~T T s b 


dAg 

35“ = 'Pc' 


- a cp + 3(cp" - a 2 cp ) + 


r + 3 c 


iaRe 

23 


(cp" - a 2 cp) + - ~ r- ■ 1 ^ 


(T + 3)' 


(29) 

(30) 

(31) 

(32) 

(33) 


(34) 
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ini nil 1 1 mi 


Sa_ 

5T * 'C + - V K * “V + 


X + cc a 


sr = <p b" + - 13 (c Pb + «pfa) + s b 


SA 3 2 r 

ST’ = ^c" + <*P” ' 3 + acp c ) + s c + iaRe 


/ . , _ aPr s 

(cp» + 2cp - r— 


2p / , \ 2 

(r + a) 


(35) 

(36) 
(3?) 


For a given point (a. Re) in the a, Re-plane, the procedure outlined should con- 
verge to an eigenvalue c. 


Criterion for Convergence to an Eigenvalue 

The procedure outlined previously for finding the eigenvalue for a given 
point in the a. Re-plane was programed for solution hy using complex arithmetic 
on the IBM 7090 located at the Lewis Research Center* The numerical integration 
was done with the Runge-Kutta method, which gives fourth-order accuracy. Con- 
vergence to an eigenvalue was established by requiring that 

| A.J 2 + 1 Agl 2 + 1 A 3 J 2 < 5X10" 6 or requiring that |Aa|, |Ab|, and |Ac| all be 
less than or equal to 5X10”^. 

The step size was externally controlled and set such that there was no sig- 
nificant change in the eigenvalue when the example was rerun with a step size 
equal to one-half of the original value. 

In order to fix the size of the solution, the eigenfunctions were normalized 
after each integration so that cp was set equal to unity at the edge of the 
boundary layer. 

The effect on the results of the choice of the value of rj for the end of 
the range of integration was examined. For air r\ was taken to be 6, and for 
water T] was taken to be 5. Variations of T] from these values for both 
cases produced little or no changes in the values of the eigenvalues. As a 
matter of fact, tj was taken to be 3 for exploratory runs that resulted in good 
approximations to the eigenvalues. 


Determination of Neutral Curve 

The results of the stability analysis are usually displayed in an a, Re- 
diagram on which the curve Cj_ = 0 is drawn. On this curve Re can be consid- 
ered a function of a. The minimum value of Re for points on this curve is 
called the minimum critical Reynolds number. The disturbances corresponding to 
points in the a, Re-plane with Reynolds numbers below the critical value will be 
damped. 
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The procedure used herein to find the critical Reynolds number consisted in 
holding Re constant and plotting the imaginary part c± of the established 
eigenvalues against a. The values of a for -which Cj_ = 0 could be read from 
such a graph establishing points on the neutral curve. For those values of Re 
below the critical Reynolds number, the plot of against a -would not give 

any c± = 0. However, after the critical Reynolds number had been bracketed, a 
was held constant, and successive linear interpolations led to the critical 
Reynolds number. Other points on the neutral curve were obtained by plotting 
o± against either a or Re and holding the other constant, depending on 
which choice proved to be more convenient. Although this procedure is straight- 
forward and easy to apply, it is an indirect way to establish points on the 
neutral curve. This disadvantage can be attributed to the choice of the complex 
number c as an eigenvalue. The advantage of using c as an eigenvalue will 
become apparent after examining the differential equations to be integrated. 

Fixing a and Re and determining the relation between the four param- 
eters a, Re, and (c r , c^) by satisfying the boundary conditions (eqs. (ll) 
to (13) ) correspond to the choice of (c r , c-j_) as the unknown eigenvalue pair and 
lead to a system of 24 first-order complex differential equations to be inte- 
grated. Since points on the neutral curve (c-j_ = 0) are of primary interest, 
another choice is to set c-j_ = 0 and then attempt to find the proper relation 
among the parameters a. Re, and c r by solving equations (17) to (19) with 
c± = 0 and Acj[ = 0. Setting c± = 0 leads to the solution of six real equa- 
tions in five real unknowns. In general, this set of equations will not be con- 
sistent unless the parameters a. Re, and c r are chosen properly. However, 
the proper value of these parameters is the information being sought. Since the 
value of Cj_ has been specified, it is no longer permissible to fix both 
a and Re, and the variation of one of these parameters has to be allowed for 
in order to provide a consistent set of six real equations in six real unknowns. 
However, the introduction of this additional unknown requires the integration 
of an additional set of six first-order complex differential equations. 

Hence, the choice of (c r , Cj_) as the unknown eigenvalue pair among the pa- 
rameters a. Re, and (c r , c-j_) leads to the fewest number of differential equa- 
tions to be integrated. This certainly justifies the choice of c as an eigen- 
value for the case when the differential equations are solved using complex 
arithmetic. If the system of differential equations is written in real form, 
the number of equations to be integrated adds up to 84. However, 72 of these 
equations involved partial derivatives with respect to the starting values and 
the eigenvalue c. The number of equations to be integrated can be reduced if 
use is made of the property that the solutions of equations (l) and (2) are 
analytic functions of r\ and depend analytically on the parameters a. Re, 
and c. This means that the information provided by the integration of the 72 
real equations given by equations (20 ) to (28) can be obtained by the integra- 
tion of only 36 real equations and the use of the Cauchy-Riemann relations, thus 
giving a total of 48 real equations to be integrated. For example, the rate of 
change of cp with respect to is given in terms of derivatives with respect 

to c r by relations of the form: Scpi/Sc^ = Scp r /Sc r and S<p r /Bcj_ = -5cpj_/Sc r . 

Insisting that Ac^_ = 0 would amount to rejecting the information furnished by 
the Cauchy-Riemann relations. 
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Another justification for the choice of (c r , c-^) as the unknown eigenvalue 
pair, thereby allowing Cj_ to take on nonzero values, is that this choice may 
lead to a solution of the eigenvalue problem, whereas setting c^ = 0 at the 
outset may eliminate any solutions. For example, the eigenvalue problem in the 
case of plane Couette flow has no solution for Cj_ = 0. 


ENERGY BALANCE OF DISTURBANCE MOTION 


After the eigenvalue problem has been solved and the eigenfunctions have 
been determined, the energy balance of the disturbance motion can be computed. 

The time rate of increase of the disturbance kinetic energy per unit of 
volume of a fluid particle that moves with the basic flow is 


JLfe (u 2 + V 2 )] = (± + u JL\[p 

DtL 2 J [hr S*/ L 2 


<V/2 v 

(u + V ) 


(38) 


An equation governing the time rate of increase of the kinetic energy can 
be obtained by the same technique as used in reference 8 for ordinary boundary 
layers, that is, by multiplying equation (B6) by u and equation (B7) by v, 
adding these equations, and using equation (B5) to simplify. This procedure 
yields 



where £ = (dv/dx) - (du/dy). 


Of course, all disturbance quantities are taken to be real in the energy 
balance equation (39). The terms in this equation are to be integrated with 
respect to y from y = 0 to y = » and with respect to x ovSr a wavelength 
A of the disturbance in order to obtain the growth of kinetic energy per unit 
of time depth and area. After integration the second and last terms on the right 
side of equation (39) vanish since u and v vanish for y = 0 and y = <» 
and have the period A with respect to x. The equation governing the time rate 
of increase of kinetic energy of the disturbance motion is 




ay 


~ ~ an 

-uv — + 

ay 


g0ut 



( 40 ) 
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The first term in the integrand of the right side of equation (40) gives the 
rate at which the basic flow, in virtue of its shearing motion, is working 
against the Reynolds stresses arising from the disturbance, the last term repre- 
sents the dissipation of the disturbance motion, and the second term represents 
the work done by the buoyancy force. 


The reference quantities introduced in appendix B are used to transform 
equation (40) to the following dimensionless form: 


JL KE = 
Dt 


Jf 

|=o *A]=o 


d| dr) 


U V pi" + _1_ _U_ _l_/5 V _ 5 u\ g 

U* U* Re AT U* Re\5I U* ~ 5^ U */ 


(41) 


where 


KE = 


// 

^1=0 J T]=0 


( aj9 rwP 

I p 


(42) 


The disturbance quantities in the integrand on the right side of equa- 
tion (41) are given by the following relations: 


u. 

IT 


~ = £e|cp' (r\) exp |ia(£ - ct)]| 
^ = &e j-iacp(r|) exp [ia(i - ct)]| 
= fe|s(rj) exp [i<x(£ - ct)]| 


t 

AT 


where fte denotes the real part of the complex quantity. 


(43) 

(44) 

(45) 


The integrals appearing in equation (4l) are evaluated through substitution 
of equations (43) to (45) followed by integration with respect to It is 

found that, for neutral disturbances, 

/ DO >*00 

d-tj (cp ' cp . - cp'cp )F" + _i_ / dri(s cp* + s.cp' ) 

lv Vi aRe / r Y r rV 

•'n 


db / d11 [ (a2cp r - <pp 2 + (« 2< pi - 


(46) 


The following designations are used to identify the integrands in equation (46): 


II 



( 47 ) 


= -(cp'cp - cp'cp )F" 
Re ^Vi 4 y r 


e B = db (S X + 8 i < Pi ) 


(48) 


S D ~ oRe 


- cpp 2 + (a 2 cp ± - cp'^) 2 ] 
For neutral disturbances the energy balance gives 


(49) 


/ OO >*00 >*00 

e Re + J o 4-1 e D - 0 


(SO) 


The satisfaction of equation (50) provides a check on the solution of the eigen- 
value problem. 


BASIC VELOCITY - AM) TEMPERATURE PROFILES FOR WATER 

The Pr = 0. 733 velocity and temperature profiles are tabulated in refer- 
ence 9. However, for the case Pr = 6.7 results were not available. In order 
to obtain the Pr = 6.7 velocity and temperature profiles, the free-convection 
boundary value problem was solved by treating it as an initial value problem as 
described in the INTRODUCTION. However, in this case there is no eigenvalue in 
the equations, rather the proper starting values have to be determined in order 
to satisfy conditions at infinity. The boundary value problem as given in refer- 
ence 9 consists of solving 

F m + 3FF" - 2(F') 2 + H = 0 (51) 

H” + 3PrFH' =0 (52) 

subject to the boundary conditions 

F(0) = F’ (0) = 0 H(0) = 1 (53) 

F' (eo) = H(«>) = 0 (54) 

The tabulation of the Pr = 6.7 velocity profile is given in table I and is 
shown in figure 1 along with the Pr = 0.733 profile. 

With the velocity and temperature profiles for both air and water available, 
the stability of these two free-convection boundary- layer flows can be investi- 
gated. 


RESULTS AND DISCUSSION 

Stability calculations were carried out for the free convection of air 
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(Pr = 0.733) and ¥ate r (Pr = 6.7) with and -without temperature fluctuations. 


Stability Results for Air Without Temperature Fluctuations 

Neglecting temperature fluctuations amounts to solving the ordinary 
boundary- layer stability problem -with the prescribed free-convection velocity 
profile inserted in the Orr-Sommerf eld equation, which is equation (l) with s* 
set equal to zero. 

The results of the present method can be compared with the calculated re- 
sults of reference 5 for the Pr = 0.733 velocity profile. Figure 2 shows the 
neutral curve drawn in the a, Re-plane calculated by the present method and also 
shows some of the neutral points calculated in reference 5. Also shown in this 
figure is the point at which a finite disturbance was observed at Re = 400 and 
a = 0.367 (ref. 10). The present method gives a minimum critical Reynolds num- 
ber of 105 at a = 0.4. This value compares well with the value of 103 calcu- 
lated in reference 5. The lowest value of Re calculated in reference 2 for 
this case was Re = 478 at a = 2.54. The values obtained from references 5, 

10, and 2, respectively, were all converted to the values quoted, which corre- 
spond to the choice of reference quantities used herein. 

In figure 3(a) the eigenfunctions cp r/ cp^, cp.jL, and cp^ are shown for the 
minimum critical Reynolds number of 105 at a = 0.4 with c r = 0.1513. Also 
shown in this figure are the locations of the critical layers. The qualitative 
agreement between these curves and the curves given in reference 5 for this case 
is good. 

The distribution of the energy transfer functions throughout the boundary 
layer is also presented in figure 3(b), where the location of the critical 
layers is also shown. As is the case for ordinary boundary layers, most of the 
energy transfer from the basic flow by the Reynolds stress takes place at a 
critical layer. Of course, for ordinary boundary- layer velocity profiles there 
is only one critical layer. For the free-convection velocity profile, the outer 
critical layer appears to be more significant than the inner critical layer. As 
can be seen from figure 3(b), the dissipation is greatest near the plate. Of 
course, for neutral oscillations the net area -under the two curves should add up 
to zero. Integration by Simpson f s rule verified this for this case and for all 
cases to be discussed subsequently. 


Stability Results for Water Without Temperature Fluctuations 

The results for this case are qualitatively the same as for the previous 
case. Figure 4 shows the neutral curve drawn in the a,Re-plane calculated by 
the present method and also shows the points at which finite disturbances were 
observed in reference 1. These experimental points were obtained from figure 8 
of reference 1. As can be seen from figure 4, all the experimental points lie 
within the region of amplification as obtained by the present method. The ex- 
perimental points represent natural finite disturbances, which were observed 
without the use of some device to introduce controlled disturbances in order to 
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provoke the onset of turbulence. The lowest value of Re calculated in refer- 
ence 1 for this case was Re = 5040 at a = 1.5. 

In figure 5(a) the eigenfunctions cp r , cp.j_, cp^., and qp£ are shown for the 
minimum critical Reynolds number of 385 at a = 0.4125 with c r = 0.0800, and 
the distribution of energy transfer throughout the boundary layer is shown in 
figure 5(b). For the case of no temperature fluctuations, the plots of the 
eigenfunctions and the distribution of energy are quite similar for air and 
water. 


Stability Results for Air With Temperature Fluctuations 

In reference 1, it is indicated that for large values of aRe the tempera- 
ture fluctuations can be neglected for the purpose of solving the eigenvalue 
problem. Of course, "large" is a relative term. For computational purposes, a 
large Reynolds number means that the asymptotic method may be applicable. In 
another sense, the statement that the Reynolds number is not too large can mean 
that the direct numerical method is applicable. Of course, the regions of appli- 
cability of these two methods could overlap. Since no numerical difficulties 
were being encountered in the stability calculations with no temperature fluctu- 
ations, it was decided that the Reynolds number was not too large. However, the 
task remained of determining whether or not the Reynolds number was large enough 
to neglect temperature fluctuations. For this reason, it was decided to examine 
the effects of temperature fluctuations. 

Figure 6 shows the neutral curve drawn in the a. Re-plane for this case. 

It is seen that the neutral curve develops a hump for small values of a and Re 
and the minimum critical Reynolds number is lowered. The calculated point 
closest to the minimum critical point is at Re = 64, a = 0.15 with 
c r = 0.2692. Also shown in this figure is the neutral curve for no temperature 
fluctuations. If a path along the neutral curve for the case with temperature 
fluctuations is followed, proceeding to lower values of a below 0.15, it is 
seen that c r increases continually until, at Re equal to about 70, the phase 
velocity of the disturbance wave is greater than the maximum velocity of the 
basic flow (see fig. l) . Consequently, there are no critical layers for values 
of a less than about 0.13. This result is impossible for ordinary boundary- 
layer flows for which the Orr-Sommerf eld equation is applicable (see ref. ll). 

It is of interest to examine the eigenfunctions and energy distribution, 
while proceeding along the neutral curve to lower values of a with c r in- 
creasing. Figures 7(a), 8(a), and 9(a) show the eigenfunctions, and fig- 
ures 7(b), 8(b), and 9(b) show the energy distributions for a = 0.45, 0.15, 
and 0.04, r esp e ct ively . 

Examination of figure 7(a) shows that the plots of the velocity eigen- 
functions with temperature fluctuations resemble the plots of the velocity eigen- 
functions for no temperature fluctuations (fig. 3(a)). Also shown in figure 7(a) 
are the temperature eigenfunctions. As the values of a become smaller, the 
plots of the velocity eigenfunctions change their shape, and the curves tend to 
oscillate less than at higher values of a (see figs. 7(a), 8(a), and 9(a)). 
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Examination of the energy distributions for a = 0.45 (fig* 7(h)) shows 
that the outer critical layer is significant in that the peak of the Reynolds 
stress-term curve is located near it, as is the case without temperature fluctu- 
ations. It can he seen that the energy distribution is qualitatively the same 
as for no temperature fluctuations, except that the temperature fluctuation term 
is reinforcing the Reynolds stress term. 

When the two critical layers are almost coincident at a ~ 0.15, it is seen 
in figure 8(h) that the critical layers have no correlation with the energy peak 
hut the buoyancy term is giving a positive contribution as in the previous case. 
Hot*, that in this case the Reynolds stress term is not adding energy to the dis- 
turbance flow, but is actually subtracting energy from it. When the critical 
layer has completely disappeared as in figure 9(b), it is seen that the buoyancy 
term is still adding energy to the disturbance flow, and again the Reynolds 
stress term is subtracting. 

From the previous results it appears that the neglect of temperature fluctu- 
ations for the purpose of solving the eigenvalue problem can be justified in the 
vicinity of a = 0.35 to 0.55, where the eigenvalues for the case without tem- 
perature fluctuations are close to the eigenvalues for this case. However, as 
can be seen from figure 7(a), the temperature fluctuations are not negligible 
even in this range of a. 

For values of a smaller than a = 0.35, the buoyancy term in the energy 
balance is significant in that It is the only term that is giving a positive 
contribution to the energy of the disturbance motion. For a = 0.15 and 0.04 
it can be seen that the Reynolds stress term is actually extracting energy. It 
appears that the Introduction of temperature fluctuations introduces a new mode 
of instability. This new mode is characterized by the buoyancy force term in 
the energy balance assuming a dominant role. The buoyancy force term is inde- 
pendent of any property of the basic flow velocity profile and Is proportional 
to the gravitational acceleration. 

The presence of buoyancy effects allows the phase velocity of the disturb- 
ance wave to be greater than the maximum velocity of the basic flow. In other 
problems of hydrodynamic stability where the energy supply to the disturbance is 
proportional to the gravitational force, the phase velocity is greater than the 
maximum velocity of the basic flow. In fact, in reference 12 it is shown, for 
small wave numbers, that the phase velocity is equal to twice the maximum veloc- 
ity for laminar flow down an inclined plane. 


Stability Results for Water With Temperature Fluctuations 

The amount of information obtained for this case is significantly less than 
for the previous case of air with temperature fluctuations, because numerical 
difficulties were encountered for large values of a and Re. Figure 10 shows 
the neutral curve drawn in the a. Re-plane. All the points on this curve repre- 
sent eigenvalues which have the property that the phase velocity of the disturb- 
ance wave is greater than the maximum velocity of the basic flow. For air 
(fig. 6) this property was obtained only for a < 0.13. Also shown in figure 10 
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is the neutral curve for water without temperature fluctuations. It was diffi- 
cult to trace the curve with temperature fluctuations for higher values of a 
than shown, because the values of the eigenfunctions at the edge of the boundary- 
layer changed markedly with each run even though the eigenvalues and starting 
values were only changing in the eighth decimal place. However, even though the 
numerical method cannot produce the eigenfunctions, the eigenvalues so obtained 
can be accepted as being reliable (see ref. 7). The point at a = 0.75 and 
Re = 385 in figure 10 is one where eigenvalues, but not the eigenfunctions, can 
be obtained. This point is on a neutral curve with phase velocity c~ = 0.1241, 
which is less than the maximum velocity in the basic flow (see fig. 1). 

The eigenfunctions for Re = 34 at a = 0.45 with c r = 0.1556 are shown 
in figure 11(a), and the energy distribution for this case is shown in fig- 
ure 11(b). It can be seen that the two cases, air and water, with temperature 
fluctuations are qualitatively the same in that the buoyancy term is providing 
most of the energy input into the disturbance motion when the phase velocity is 
greater than the maximum velocity in the boundary layer. 


CONCLUDING REMARKS 

The direct numerical method gives a minimum critical Reynolds number lower 
than the Reynolds number at which finite disturbances were observed experimen- 
tally, whereas calculations based on asymptotic techniques yield a minimum 
critical Reynolds number higher than the Reynolds number at which finite dis- 
turbances were observed. However, the direct method is not universally appli- 
cable to all problems for all ranges of the parameters, especially when the 
Reynolds number is large. Generally, the method should give results in the 
lower left corner of the a, Re-plane where the minimum critical Reynolds number 
is usually found. 

By referring to the results of experiments with natural disturbances, that 
is, experiments conducted without the use of controlled disturbances to provoke 
the onset of turbulence, it can be concluded that the minimum critical Reynolds 
number calculated herein for the two cases without temperature fluctuations pro- 
vides a lower bound to the Reynolds number at which finite disturbances were ob- 
served. However, inclusion of temperature fluctuations indicated that this mini- 
mum critical Reynolds number is not the least lower bound. Rather, there is a 
lower minimum critical Reynolds number at small values of wave number a. This 
instability, for which the phase velocity of the disturbance wave is greater than 
the maximum velocity of the basic flow in the boundary layer, has not been re- 
ported in the accounts of natural transition experiments. 

The existence of temperature fluctuations provides a new mechanism of 
energy transfer to the kinetic energy of the disturbance motion. This amount of 
energy is a significant contribution to the energy balance in that, for all 
cases calculated, it was always positive, that is, destabilizing. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, September 19, 1963 
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.APPEBDIX! A 
SYMBOLS 

a cp'" (0) 

b s’ (0) 

c phase velocity 

c specific heat at constant pressure 

ir 

e B buo;yancy 

e-Q dissipation 

e^ e Reynolds stress 

F stream function, dimensionless 

Fr Froude number 

Gr^ Grashof number based on x, gp* ATx 3 /v^ 
g gravitational acceleration 

H temperature dimensionless 

i imaginary unit 

KE defined in eq_. (42) 

k coefficient of heat conductivity 

P pressure of basic flow 

Pr Prandtl number 

p pressure 

Re Reynolds number 

real part 

s temperature amplitude function, dimensionless 

T temperature, basic flow 

- T 

w 00 
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wall temperature 
^ ambient temperature 

t temperature 

U reloelty parallel to plate,, basic flow 

U* reference velocity parallel to plate 

u velocity parallel to plate 

V velocity normal to plate, basic flow 

v velocity normal to plate 

x distance from leading edge of plate 

y normal distance from plate 

a wave number 

0 * coefficient of volumetric expansion 

5 reference length, -y/2x/ 

T) normal distance from plate, dimensionless 

tj edge of boundary layer 

7 \ wavelength, 2rt/a, dimensionless 

p coefficient of viscosity 

v kinematic viscosity 

| distance from leading edge of plate, dimensionless 

p density 

t time, dimensionless 

9 stream function amplitude, dimensionless 

i|c stream function, dimensionless 

Subscripts: 

a, b, c denotes differentiation with respect to quantity 
i refers to imaginary part 
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max 


maximum 


r refers to real part 

Superscripts: 

~ disturbance quantity 

— dimensional quantity 

1 differentiation with respect to r) 
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APPENDIX B 


DERIVATION OP STABILITY EQUATIONS 

The governing equations for a parallel basic flow plus a disturbance as 
given in reference 1 are as follows: 


5u 

St 


u 


5u 

57 

5v 

St 


5u , 5v _ n 
57 + 57 

(Bl) 

= - i ^ + v^u + g p*(t - tj 

(B2) 

^Z+v^Z=-i^P + v^v 

OX dy p dy 

(B3) 

■ u ^.+ v^i = -n.v 2 t 

dx dy pc 

(B4) 


where 


v 2 = ^L + lL 

5x 2 5y 2 

The "basic steady-state flow is the f ree-convect ion flow about a vertical heated 
plate. In the derivation of the governing equations the variation of the density 
has been neglected except in the gravity term where the density was taken to de- 
pend linearly on the temperature. In equations (Bl) to (B4) the density, where 
it appears, is taken to be a constant. In accordance with the parallel-flow 
assumption, the basic flow can be described by U = U(y), V = 0, P = P(y), and 
T = T(y). Superimposed upon the basic flow is a two-dimensional disturbance. 

The disturbance equations are obtained by substituting u = U + u, v = v, 
p = P + p, and t = T + t in the previous equations, subtracting out the basic 
flow, and neglecting products of disturbance quantities and their derivatives. 

The disturbance equations are 



5u . 5v _ o 
57 57 

(B5) 

St 

u 5u + ~ 5u = _ l 5p + vV 2 ~ + g3 ^ 

dx dy p dx 

(B6 ) 


*1 + U = - i p + vV 2 v 

57 ox p dy 

(B7) 


a* + u + 7 dT = k V 2 ! 
5~ ^ dy pc p 

(B8) 
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The disturbance velocities can he obtained from a stream function 

t(x,y,T) = cp(y) exp |ia(x - ct)J (B9) 

from -which 

U = ^ = <p* (y) exp [ia(x - cr)] (BIO) 

v = - = - iacp(y) exp [ia(x - cT)J (Bll) 

As can he seen from equation (339), the disturbance is taken to he periodic 
in the distance x along the plate. The positive quantity a is the wave num- 
ber of a disturbance wave, and c* r , the real part of c, is the velocity of 
propagation of the wave. The imaginary part of c* will determine whether the 
disturbance will grow > 0) or decay (c*-^ < 0) in time. The temperature dis- 

turbance is also periodic in the distance x along the plate and may be ex- 
pressed as 

t(x,y, t) = a (y) exp [ia(x - ct)] (B12) 

It is convenient to represent the disturbance quantities in complex form in order 
to satisfy the phase relations imposed by equations (335) to (338). However, 
physical significance is to be attached only to the real part of disturbance 
quantities. 


Dimensionless variables are introduced by choosing reference quantities in 
conformity with those used in calculating the basic flow as follows: 

Length: 


Velocity: 


Temperature: 


where 


yfo . = & 


2v(Gr x ) 1 / 2 

x 


= U* 


T - T - AT 1 
v °° 


Gr = gP* 
x 2 

v 
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The following dimensionless -variables are used to transform equations (B5) to 
(B8) into the dimensionless form: 


x = £5 

(B13) 

I! 

-rJ 

o> 

(B14) 

_ t5 
T _ U* 

(B15) 

U = F'U* 

(B16) 

- T = HZ!T 

CO 

(B17) 

9 = 95U* 

(B18) 

s’ = s AT 

(B19) 

— TT* 

G = CU 

(B20) 

— a 

a = 8 

(B21) 


In terms of the new dimensionless variables, the stability equations are obtained 
by substitution of equations (BIO) and (Bll) into equations (B5) to (B8); the 
result is 

(F' - e)(cp" - a 2 cp) - F'"cp — (cp"" - 2a 2 cp" + a 4 cp) - j-.§* s' (B22) 

aRe aFr 

(F 1 - c ) s - cpH- = - — i— (s" - a 2 s) (B23) 

aRePr 

where the Reynolds number Re is defined as Re = SU*/v = 2-/2(Gr x )‘*" / ^, the 
Froude number is Fr = (U*) 2 /gS = 4v 2 Gr x /g8x 2 , and the Prandtl number is 
Pr = qCp/k. Equations (B22) and (B23) become identical to equations (l) and (2), 
respectively, when the substitution Fr/|3* AT = Re is made. 

Equations (B22) and (B23) differ from the ones given in reference 2 in that 
the term s' in equation (B22) is multiplied by F 2 / (F^ ny ) 5 . This difference 
is due to the choice of reference quantities in reference 2, which differ from 
those of reference 1 used hereinafter. The reference quantities of reference 2 
are as follows: 


Length: 


F °° V2x 

F max (Gr ) 1//4 
x 
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Velocity: 


F f 

max 




APPENDIX c 


CONDITIONS TO BE SATISFIED BY EXPONENTIALLY 
DECAYING SOLUTIONS 

The technique used in reference 4 is used to eliminate the arbitrary con- 
stants in equations (7) and (8) in order to obtain a set of linear and homogene- 
ous relations in the dependent variables of equations (l) and (2) and their de- 
rivatives, which are to be satisfied at the edge of the boundary layer. 

In order to satisfy the boundary conditions (eq. (4)), it is necessary that 
ci = c 3 = c 5 = 0. Further, in order that the solutions (7) and (8) agree with 
the numerical solution of equations (l) and (2), which is obtained by numerical 
integration when q = q , it is necessary that 

9 = Cq exp (aq g ) + c g exp (-aq g ) + exp (Pq g ) + c 4 exp (-0q e ) 

+ c,- exp (m e ) + c g exp (-rq e ) ( C1 ) 
9' = c-ja exp (aq ) - c 2 <x exp (-aq e ) + exp (0q g ) - c 4 P exp (-0q 3 ) 

+ c 5 r exp (rq e ) - e g r exp (-rn 3 ) (02) 

cp" = CqO 2 exp (aq e ) + c 2 a 2 exp (-aq g ) + c 3 P 2 exp (Pq e ) + c 4 P 2 exp (-0q e ) 

n O 

+ c 5 r exp (rq e ) + c g r exp (-rq e ) (C3) 
cp'” = Cl a 3 exp (aq e ) - Cgtt 3 exp (-aq e ) + c 3 (B 3 exp (0q e ) - c 4 p 3 exp (-pq g ) 

+ c r 3 exp (rn ) - e r 3 exp (-m ) (C4) 

O So 0 

(r 2 - a 2 )(r 2 - P 2 ) exp (rq ) (r 2 - a 2 Hr 2 “ P 2 ) exp ( -yq ) 

s = - c_ + c c , , e - (C5) 

D y O y 

s’ = -c 5 (r 2 - a 2 )(r 2 - P 2 ) exp (rq e ) - c g (r 2 - ct 2 )(y 2 - p 2 ) exp (-rq g ) (C6) 

where the left sides of equations (Cl) to (C6) are obtained from the numerical 
integration of equations (l) and (2) from q = 0 to q = q e . Solution of equa- 
tions (Cl) to (C6) gives the c's. Since Cq = c 3 = c 3 = 0, the three determi- 
nants that give Cq, c 3 , and c g must vanish. Hence, 
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9 


9' 


s' 


1 1 


cp' -a 3 


I cp" a 2 


r 3 


0 0 


1 

-3 

2 


-a 3 3 5 - 3 3 


0 0 


r 

2 \ / 2 


/ 2 2 \ / 2 0 2 \ 

(r - a )(r - P ) 


-r 

r 2 

*2 

-IT 


_(r* - a a )(r z - 3 g ) (r 2 - a 2 )( r 2 - 


3 2 ) 


r 


■(r 2 - a 2 )( r 2 


3 2 ) 


9 


a -a cp 1 


a 


a 


a 2 cp" 


1 

-3 

,2 


-a 3 cp"' -3 3 


1 

r 

r c 

*2 

r" 
. 2\ /„2 


_ (r - ar) (y~ - 3^) 


r 


0 s' 0 - (r* - a 2 )(y 2 - 3 2 ) - (t 2 - cl 2 )(y‘ 2 - 3 2 ) 


1 

-r 

r 2 

-T 5 

(r 2 - a 2 )(r 2 


- 3 2 ) 


2w_2 


111 

a -a 3 


1 

■3 


2 2 2 0 2 
a a 3 3 


9 

cp' 

cp" 


5 3 „3 3 

a -a 3 -3 9 


0 0 


1 

-r 

2 

r 

-2 

-r" 


/ 2 2 w 2 

(r - a ) (r 


3 2 ) 


= o 


Equations (C7) to (C9) reduce to 

(r 2 - a 2 )[cp'" + a9" - 3 2 ( 9 ' + 09)] + r 2 s + as' = 0 
(r 2 - 3 2 ) Op'" + 39" - cg 2 ( 9 ' + 39 )] + r 2 s + 3s' = o 


= 0 (C 7 ) 


= 0 (G8) 


(C9) 


(CIO) 

(Oil) 
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s' + ys = 0 


(C12) 

Equations (CIO) to (C12) are the conditions that the numerical solution must sat- 
isfy at the edge of the boundary layer at r| = T] e . It is convenient to eliminate 
s' from equations (CIO) and (Cll) by means of equation (C12). This leads to the 
somewhat simpler conditions that will be employed instead of equations (Cll) 
and ( CIO ) , respectively, namely, 

cp'" - a 2 cp' + 0(q>" - a 2 cp) +■ 8 = 0 (C13) 

9 "' + acp" - |3 2 (cp' - acp) + a = 0 (C14) 

I T CL 

It is clear that, if equations (C12) to (C14) are satisfied, then so are equa- 
tions (CIO) and (Cll). 
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TABLE X . - FUNCTIONS F AND H AND DERIVATIVES 
FOR PRANDTL NUMBER OF 6.7 


T ) 


0 . 

.1250 

.2500 

.3750 

.5000 

.6250 

.7500 

.8750 

1.0000 

1.1250 

1.2500 

1.3750 

1.5000 

1.6250 

1.7500 

1.8750 

2.0000 

2.1250 

2.2500 

2.3750 

2.5000 

2.6250 

2.7500 

2.8750 

3.0000 

3.1250 

3.2500 

3.3750 

3.5000 

3.6250 

3.7500 

3.8750 

4.0000 

4.1250 

4.2500 

4.3750 

4.5000 

4.6250 

4.7500 

4.8750 

5 . 0000 

5.1250 

5.2500 

5.3750 

5.5000 

5.6250 

5.7500 

5.8750 

6.0000 

6.1250 

6.2500 

6.3750 

6.5000 

6.6250 

6.7500 

6.8750 

7.0000 

7.1250 

7.2500 

7.3750 

7.5000 

7.6250 

7.7500 

7.8750 

8.0000 

8.1250 

8.2500 

8.3750 

8.5000 


F 


0 . 

.0032 

.0118 

.0241 

.0388 

.0548 

.0715 

.0880 

.1041 

.1194 

.1338 

.1472 

.1595 

.1707 

.1810 

.1904 

.1989 

.2066 

.2136 

.2200 

.2257 

.2309 

.2355 

.2397 

.2436 

.2470 

.2501 

.2529 

.2554 

.2577 

.2597 

.2616 

.2633 

.2648 

.2661 

.2673 

.2684 

.2694 

.2703 

.2711 

.2718 

.2725 

.2731 

.2736 

.2741 

.2745 

.2749 

.2752 

.2755 

.2758 

.2760 

.2763 

.2765 

.2766 

.2768 

.2769 

.2771 

.2772 

.2773 

.2774 

.2775 

.2775 

.2776 

.2776 

.2777 

.2777 

.2778 

.2778 

.2778 


F » 

0 . 

.0494 

.0852 

.1095 

.1244 

.1318 

.1335 

.1311 

.1259 

.1189 

.1110 

.1026 

.0943 

.0862 

.0786 

.0714 

.0648 

.0587 

.0532 

.0481 

.0435 

.0393 

.0355 

.0321 

.0289 

.0261 

.0236 

.0213 

.0192 

.0173 

.0156 

.0140 

.0127 

.0114 

.0103 

.0093 

.0083 

.0075 

.0068 

.0061 

.0055 
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Figure 6. - Neutral curve for air with and without temperature fluctuations. 
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